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abstract 



We investigate the contributions of neutral Higgs bosons to nonleptonic transition b -^ 
^ ■ qss, q = d, s under the supersymmetric context. Their effects to decay width and CP 

cn . violation in corresponding exclusive decays are explored. The anomalous dimension matrices 



of the operators which have to be incorporated to include the contributions of neutral Higgs 



^>0 ', bosons are given. We find that when tan/? is large (say, 50) and neutral Higgs bosons are not 

too heavy (say, 100 GeV), contributions of neutral Higgs penguin can dominate electroweak 
Q\ • penguin contributions, and for some processes, they can greatly modify both decay width 

r-| ■ and CP asymmetry. 
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The systematic study of the nonleptonic decays of B mesons has been led long before the 
construction of B factories. The CLEO report [1] on the measurement of the charmless two- 
body hadronic decays revived the theoretical interest in this subject. Among the reported 
exclusive decay modes, B -^ Krj' received much more attention (for a review, see [2]), 
and according to [3, 4], four-quark-operator contribution is capable of explaining it under 
the theoretical framework of SM. For the inclusive decay B -^ Xgi]', it seems that new 
mechanism or new physics should be introduced in order to explain the unexpected large 
branching ratio [5]. 

The effective Hamiltonian method is the fundamental tool to study the charmless non- 
leptonic decays of B mesons. Starting from a full theory, after integrating out the heavy 
freedom degree at a certain energy scale, the effective Hamiltonian is obtained. For the 
standard model (SM), the heavy freedom degrees integrated out at /iq = m^ are top quark, 
Z boson and W boson. For the SUSY models, in order to be simple, the extra heavy freedom 
degrees integrated out include all superpartners, though there are superpartners, such as the 
lightest neutralino, which may be lighter than W boson, and there are superpartners, such 
as the super-quarks of the first two generations, which may be much heavier than W boson. 
At this integrating out scale /iq = rriw, super-partners can not only contribute to the effective 
Wilson coefficients, but also introduce new operators which will be shown below in the case 
of including contributions of neutral Higgs bosons (NHB). 

The works on the next-to-leading-order QCD-improved effective Hamiltonian derived from 
SM were reviewed in [6], where the problem on the scale independence were studied. To 
present days, our knowledge of the dynamics of the hadronic transitions are limited, and 
the approach employed in evaluating hadronic matrix elements of operators is factorization. 
There are two last papers which represent the state of art of the systematic study of charmless 
hadronic two-body B decays. Ref [3] concentrates on the experimental test on the naive 
factorization ansatz, having classified 76 decay modes into five classes and analyzed the 
dependence of these decay modes on the weak mixing matrix elements, form factors, decay 
constants, QCD scale, and the effective number of color. Ref [4] shows the experimental 
data of p°7r^ and (pK^ cannot be accommodated simultaneously by treating Nl^f{LL) = 
Nff{LR) and, by assuming that Nff{LL) ~ 2 and Nff{LR) ~ 5, finds the existed 
experimental data can be explained well. 

Nonleptonic decays have also been investigated in SUSY models [7, 8]. Several papers deal 
with the possible effects of SUSY to hadronic B decays in the SUSY models with R parity 
violation [9]. The other topic related with the nonleptonic B decays of SUSY models is CP 
violation, there are a vast amount of literatures [10, 11]. 

It is well known that for charmless B hadronic decays with Ai? = 1 in which tree level 
contributions are suppressed by CKM and color symmetry, penguin contributions, among 
which gluonic ones play major parts, always dominate. In this letter we investigate the 
total contributions of penguins in MSSM, including NHB penquins in addition to QCD 
penguins and electrweak penguins. The motivation of this paper is based upon the fact 
found by us in [12], where it is found that in SUSY models, when tan/5 is large, NHB 
are not heavy, and the mass splittings of stops and charginos are large, contributions from 
exchanging NHB can greatly modified the branching ratio and backward-forward asymmetry 
of i? — ^ Xgl^l" .,1 = fi,T. It is a natural further step to explore their contributions to 
nonleptonic B decays. Another motivation for our work is the striking large branching 



ratio for the exclusive decay of -B ^ Krj' [13], which has caused much theoretical efforts to 
assess the theoretical prediction of SM and the validity of factorization approach [14, 3, 4]. 
Our results show that contributions of NHB can be larger than the sum of the electroweak 
contributions from exchanging 7 and Z bosons, and for some decay modes, they can modify 
both the decay widths and CP asymmetries significantly. 

We shall concentrate on the decay modes in which the quark level process is 6 — > qss^ q = 
d, s in order to see the contributions from NHB. The transition b -^ qss is purely induced 
by FCNC penguin diagrams, which can be divided into four classes, l)those of exchanging 
gluon, 2) of exchanging photon, 3) of exchanging Z, and 4) of exchanging Higgs neutral 
bosons. Normally, contributions from the first class are expected to be large, while those of 
the second and the third are small because of aem/c(s and a2/as respectively. In the SM, 
since suppressed by the ratio of masses of quarks and masses of NHB, those of the fourth 
are also negligible , so only the contributions of the first are considered. But many papers 
[15] (for a review, see [16]) have shown that contributions from the second and the third 
classes can play an important role in both the width decays and direct CP violation in some 
decay modes, and they can even play the major part in some processes, Bg — > K~^K~, for 
instance. Furthermore, in SUSY context, when tan/5 is large and the masses of NHB are not 
too heavy, contributions from exchanging neutral Higgs bosons can be considerably large 
and therefore can greatly modify theoretical predictions for some processes. So in the letter, 
we consider the case of large tan/3 and take into account all contributions of the four classes. 

The b transition b ^ q,q = d, s; can be described by the AB = 1 effective Hamiltonian 
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The first two operators are induced by tree diagrams, Oj=3,...,6 by gluonic penguins, Oj=7,...,io 
by 7 and Z penguins, Oj=ii,...,i4 and Oi=i7,i8 by neutral Higgs penguins, and Oi9,20 by chro- 
mo magnetic and magnetic penguins^. 



■^In order to make the minimal revisions we use the labels of operators as those in the original version of 
the paper. Comparing with those in refs. [26, 27, 28, 29], On, Ois, O19, O20 correspond to O15, Oig, Osg, O77 
respectively. 



The boundary conditions of the first ten Wilson coefficients and those of Oi9,2o for the 
solutions of the renormalization group equations are given in [6], here we provide those of 
Oii,...,i4,i7,i8 in the large tan/? case. 
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The conventions can be found in [12]. To reach these, we have neglected some terms 
contributing minor. It is interesting to mention that the terms proportional to tan^/5 only 
come from the one Feynman diagram where firstly virtual NHB are emitted from b-quark 
line and decay into strange quark pair, then followed with flavor change caused by loops 
exchanging charginos and stops. While for virtual NHB decaying to down quark pair, because 
of the smallness of nid/ms the contributions are negligibly small. To up-type quark pair, 
SUSY contributions at most be proportional to tan^/5 and consequently also not important 
compared to those for decaying to strange quark pair. This is the reason why we only pay 
attention to the transition h — > qss. 

Coii 13 are proportional to x^^^logx^^^^,. Therefore, contributions of SUSY depend on not 
only masses difference between stops, but also those between charginos. In order to escape 
the GIM cancellation and to increase the contributions of SUSY, both large mass splittings 
of stops and charginos and large chargino masses are needed. In contrast with the process 
6 — s> 57 and the mixing of B^ system, the smaller the mass of the lighter charginos is, the 
larger the contributions of SUSY. 

For the first ten Wilson coefficients, we use next-leading-order QCD corrected renormaliza- 
tion group equations [6]"^, while for the last six ones (i.e., Oii,...,i4,i7,i8), we use leading-order 
ones since we carried out calculations of the relevant anomalous dimension matrices only at 



''We do not include the mixing of On. 12 onto O3....,io which has been given in ref. [29]. We also do not 
include the mixing of Oi3,i4,i7,i8 onto Oi9^20 which is given first in ref. [30] and verified in ref. [29]. 



one-loop^. The anomalous dimension matrices of the last six operators (i.e., Oii,...,i4,i7,i8) 
can be divided into two distangled groups 
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Our calculation shows that Coi2,-,i4,i7,i8i^^b) are much smaller than Coni'mb), so it is ap- 
propriate to neglect them in our numerical analysis. 

As given in [6], the effective coefficients Oj are dependent of N^^^. In this note, we do not 
examine the dependence and limit ourself to a typical value of N^^^ , 3, in order to see the 
contributions of NHBs. For N^^^ = 3 and transition b ^ s, the first ten effective coefficients 
given from Wilson coefficients in SM are [3] 
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while for transition b ^ d, they are 
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These complex figures incorporate contributions of both the strong phases and weak phases. 
For the last 8 coefficients, 10^^ should be multiplied, and the numbers in the brackets refer 
to the conjugate processes b -^ sss and b -^ dss. 

The input parameters used in our numerical calculations are as follows. 

Wolfenstein parameters [17] for the CKM matrix are given as A = 0.81, A = 0.2205±0.0018 
[18], p = 0.12 and rj = 0.34 [19]. The scale for the running masses of quarks is set to 2.5 
Gev, and the corresponding masses are given in Table I in [3]. We use the BSW approach 
to evaluate the form factors [20], and the form factors at zero momentum can be found in 
the Table II in [3]. For the reason that the form factors are not sensitive to the pole masses 

'"'Wc would like to thank Dr. CD. Lii for his collaboration on the calculations. A complete analysis of 
mixing of the operators will be published else where. 



since only small extrapolations from Q^ = are involved in the decays B -^ hih2, we set 
the pole mass to 5.4, which is given in Table III in [3], in calculating the form factors . 
Decay constants are given in the Table V in [3]. The mixing angles of the flavor SU(3)-octet 
and singlet components are set to 6g = —21.2° and 6q = —9.2° [21]. For more detailed 
information on conventions, please refer [3]. 

We calculate Cqh and find that it can reach 0.05. So neutral Higgs bosons can dominate the 
contributions of the electro-weak penguins for some processes. In calculations the relevant 
ten supersymmetric parameters are taken as 
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300 GeV, niAo = 80 GeV, m^ = 500 GeV, 

= 300 GeV, At = 300GeV, tan(3 = 50 (9) 

= vr ^At = 0.3 

which satisfy the constraint on the SUSY phases arising from EDMs through the Barr-Zee 
mechanism [22, 23] as well as the constraint from 6 ^ 57 [12]. 

we get GoiiiiTT'b) = —0.051 — 0.027i, and consequently ali = —0.051 — 0.027i and all = 
—0.017 — 0.009i. We note that the parameters (9) satisfy the constraint from b — * 57. For 
the first ten Wilson coefficients, we find that the contributions of SUSY with the above 
parameters can be safely neglected, because they only modify the Wilson coefificients in few 
percent in most of part of the parameter space. So we shall still use eqs. (7) and (8) for 
the sake of simplicity. The complex phase of Cqh is originate from the SUSY CP phases, if 
tpAt = 0, the corresponding a^f = —0.055, no imaginary part appears. 

To evaluate the relevant transition matrix elements, we will use the naive factorization 
assumption. In the spectator model, several formula on decay amplitudes are listed below 
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and the superscript "efF' of a^ has been omitted for the sake of simphcity. 

Compared with formula given in [3], we have taken into account the contributions of 
NHB's by adding an, au and 013 in these amphtudes of decays. We have not included 
W^-exchange, W^-annihilation and spacelike penguin contributions because the W-exchange 
and W-annihilation contributions are negligible due to helicity suppression (as well as color 
suppression for W-exchange) and we do not have a reliable method to estimate the spacelike 
penguins[4]. We have used the equations of motion and Fierz rearrangement to calculate 
the contributions of scalar operators. For other formula, it is easy to derive when compared 
with the tables given in the appendices in [3, 4]. Our numerical analysis about the effects of 
NHB are based on these formula. 

For charged B^, Aqp is defined as 

^^^ " r(5-^/-) + r(i?+^/+) ^^^> 

and it reflects the magnitude of the direct CP violation. For the cases in which both B^ and 
B^ decay to the same final states, Acp is defined as 

^"" - (1 + lend +4) ^''^ 

where e = -34, a quantity reflecting effects from both the mixing and the direct CP violations. 
We omit the analysis on the CP violation for decays, B^ -^ VV, and B^ -^ K*^K^, K*^K^, 
in order to simplify the calculations, xb is set to 0.723 as given in [24]. We do not introduce 
new CP origin in the mixing of the neutral B meson. 
Our result is shown in Table 1. A number of observations are in order. 

1. Compared to the SM, the decay widths for most modes corresponding to the transition 
b -^ qss increase. For example, for B^ -^ K^rj', contributions of NHB increase the branching 
ratio up to 25% and for B^ -^ K^rj, they increase up to 80%. For B -^ K{K*)(f) {B^ -^ 
Ti^Tj') because of the cancellations between 012(011) and 04(05), the decay widths decrease. 
We see that the contributions of NHB do modify the branching ratios for all charmless 
decay modes whose amplitudes are governed by penguins for transition b — * qss, though 
for most nonleptonic decays gluon penguins dominate. Compared with the semileptonic 
decays b -^ s/+/~, where branching ratios and forward-backward asymmetry can be greatly 
modified [12] by the effects of NHB, for most charmless nonleptonic decay modes their effects 
are soften by the large values of the coefficients Oi, 04, Og. 

2. For those decay modes where electroweak penguins play an important role, such as 
B -^ K*r], B° -^ K*^{K^)ri, 5° -^ Ti^rj^B^ -^ vivlviv'), and 5° -^ cu^r]', the contributions 
of NHB can enhance the branching ratios by several handreds percent. 

3. For the modes where there are subtle cancellations between contributions in SM, such 
as i? -^ K*^K, B^ -^ K*^ri', and B^ -^ p^r]', the contributions of NHB can enhance 
the branching ratios by a factor of a few tens and greatly modify the corresponding CP 
asymmetry. 

In the effective Hamiltonian method to study hadronic decays the largest theoretical un- 
certainty arises from the calculations of hadronic matrix elements of local operators. For 
energetic two-body decays of B mesons factorization should be a good approximation be- 
cause of color transparency [25] and non- factorization effects may be parameterized by N^^^ . 
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In the factorization framework, in addition to the uncertainty due to non-factorization ef- 
fects, the other uncertainties come from the input parameters, such as form factors, decay 
constants, quark masses, CKM matrix elements, the scale at which Wilcon coefficients is 
calculated, etc. The different choices of input parameters can lead to an uncertainty of 
from a few percent to several handreds percent (dependent on modes), as can be seen by 
comparing the results for the same iV^-^-^ between ref . [3] and ref . [4] . The uncertainty due to 
different N^^-l" can also reach several handreds percent. Therefore, it is difficult to disentagle 
the non-standard physics for most of modes in Table 1. However, by comparing our result 
and the results given in [3] and [4], it can be drawn that 

A. compared to SM, a significantly large enhancement of the branching retio for B -^ K*r]' 
in MSSM is definitely existed; and 

B. the branching ratio for B -^ p^rj' {B^ -^ K*^K^) is enhanced by a factor of 40 (25) com- 
pared to SM for the same N^^^ and larger than all values corresponding to different N^^^ in 
SM given in [3, 4]. Thus, it could be possible to disentagle MSSM from SM by measuring 
the branching ratio if tan/? is large and NHBs are not heavy. 

In summary, we have analyzed the contributions of NHB to charmless nonleptonic two- 
body B decays. Our analysis shows that at some regions of the parameter space in MSSM, 
NHBs can dominate the contributions of electroweak penguins and significantly modify both 
decay widths and CP asymmetry of the decay modes which are due to the transition h -^ qss. 
Our analysis has not include CP phase of SUSY origin in mixing. After taking into the 
complex contributions of SUSY to the mixing, CP asymmetry of most of neutral B meson 
decay modes could be modified. 

/,From the above analysis, we know that the effects of NHB, for most charmless nonleptonic 
decay modes, are soften by the large values of the coefficients ai, 04, ag. But their effects can 
show up when there are subtle cancellations. Therefore, although there are huge theoretical 
uncertainties, the enhancements for some modes, such as i? — * p^rj' and B^ -^ K*^K^, 
exceed the theoretical uncertainties significantly so that it would be possible to search new 
physics effects by measuring these decay modes. 
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Table 1: Branching ratios in unit of 10 ^,and CP asymmetry in unit of %. 



Nc = ?) Branching Ratio CP Asymmetry 

Decay Modes SM SUSY SM SUSY 



B- - 


-^ K 1] 


1.8 


3.27 


-6.1 


-0.16 


B' - 


-^ K'T] 


25.2 


31.9 


4.5 


5.58 


B- - 


^ K'(j) 


13.8 


8.11 


2.0 


-4.8 


B- - 


-^ K*'ri 


1.50 


6.64 


5.4 


3.3 


B- - 


-^ K*~i 


0.89 


6.71 


23.7 


-5.4 


B- - 


-^K*~(t) 


6.95 


4.43 


2.0 


-3.7 


B- - 


^K-K' 


0.67 


0.82 


-12.3 


-9.8 


B - 


-^ TT^r] 


2.10 


2.14 


-14.0 


-13.3 


B- - 


-^ TT^T]' 


2.16 


2.56 


-13.0 


-11.7 


B - 


-^ K^K*~ 


0.46 


0.69 


-14.5 


-9.2 


B - 


-^ K*^K- 


0.0004 


0.01 


-46.7 


6.2 


B- - 


^ P'V 


6.16 


6.89 


-3.4 


-3.1 


B- - 


^ p-T]' 


7.41 


6.99 


-11.0 


-11.3 


B- - 


^ X*-K*^ 


5.73 


3.77 


-14.5 


-22.1 


B'- 


■^ K'^T] 


1.6 


3.1 


27.8 


37.5 


B^- 


^kS 


24.1 


30.5 


29.4 


33.5 


B^- 


-^ ^K^ 


13.3 


7.8 


29.3 


15.1 


B^- 


^ K*^r] 


1.43 


6.14 


30.0 


41.4 


B^- 


-^ K*''r]' 


0.43 


6.4 


24.6 


37.2 


B^- 


-^K*^(t) 


6.69 


4.27 


— 


— 



^0 ^ ^0^0 0.65 0.79 33.2 35.6 

^0 ^ ^0^ 0.32 0.71 31.3 38.8 

B» ^ ttV 0.12 0.019 26.8 -70.3 

B°^r]r]' 0.039 0.11 24.0 39.0 

B^ -^ r]r] 0.13 0.35 24.0 39.0 

^0 _^ ^Y^ 0.024 0.08 -4.3 -27.2 

50 ^ K'^K*o 0.44 0.66 — — 

B° -^ K*^K° 0.0004 0.01 — — 

B° -^ pr] 0.015 0.0076 34.2 -34.7 

B^ -^ pi 0.0028 0.11 63.1 27.3 

50 -^ ujr] 0.046 0.048 -10.5 -56.3 

B^^UT]' 0.019 0.11 -57.9 36.6 

^0 _^ j^*oj^*o 5_52 3.63 — — 
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